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Shell-model calculations for neutron-rich carbon isotopes with a chiral 
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We have studied neutron-rich carbon isotopes in terms of the shell model employing a realistic 
effective hamiltonian derived from the chiral N'^LOW nucleon-nucleon potential. The single-particle 
energies and effective two-body interaction have been both determined within the framework of the 
time-dependent degenerate linked-diagram perturbation theory. The calculated results are in very 
good agreement with the available experimental data, providing a sound description of this isotopic 
chain toward the neutron dripline. The correct location of the dripline is reproduced. 

PACS numbers: 21.60.Cs, 23.20.Lv, 27.20.-|-n, 27.30.+t 



The advances in the production of exotic beams are 
making possible the exploration of isotopes at and close 
to the neutron dripline. This is a major goal of nuclear 
physics, since these studies may reveal exotic features, 
such as shell modifications, that go beyond our present 
understanding of the structure of nuclei. In a recent Let- 
ter [l[ Tanaka et al. have reported on the observation 
of a large reaction cross section in the dripline nucleus 
^^C. The structure of ^^C is quite interesting, because it 
can be considered the heaviest Borromean nucleus ever 
observed [2|- This nucleus is weakly bound with a two- 
neutron separation energy S2n evaluated to be 420 ± 940 
keV Jj], and may be seen as composed of three parts: two 
neutrons plus ^°C. These three pieces must all be present 
in order to obtain a bound ^^C because of the particle 
instability of ^^C. This is the structure of a Borromean 
nucleus, which has only one bound state but, considered 
as three-body system, admits no bound states in the bi- 
nary subsystems Q. 

The scope of the present work is to perform a theo- 
retical study of the properties of heavy carbon isotopes 
approaching the neutron dripline. Our framework is the 
nuclear shell model, with both single-particle energies 
and residual two-body interaction derived directly from 
a realistic free nucleon-nucleon [NN) potential. This 
fully microscopic approach to shell-model calculations, 
where no adjustable parameter is introduced, has al- 
ready proved to be successful in reproducing the neutron 
dripline for the oxygen isotopes [J] as well as the spec- 
troscopicproperties of the p-shell nuclei @ and TV = 82 
isotones [a|- 

This is the first time that a realistic shell-model calcu- 
lation is performed in this region. In fact, most of shell- 
model calculations to date have been performed employ- 
ing empirical single-particle energies and two-body ma- 
trix elements adjusted to reproduce selected experimen- 
tal data using the full psd model space (see for example 
[3, i_S] and references therein). 

We now give a short description of our calculations. 
We have carried out shell-model calculations, using the 



Oslo shell-model code [9i|, for the heavy carbon isotopes 
in terms of valence neutrons occupying the three single- 
particle levels 0(^5/2, 0^3/2, and lsi/2, with ^^C consid- 
ered as an inert core. This choice is mainly motivated by 
the observed large energy gap {^ 6 MeV) between the 
ground and first excited state. It is worth mentioning 
that in [T^] it is suggested that the anomalous suppres- 
sion of the observed B{E2;2^ -^ 0^) transition rates 
in i4-i8(^ indicates a possible proton-shell closure in the 
neutron-rich carbon nuclei. 

The shell-model effective hamiltonian HcS has been de- 
rived within the framework of the time-dependent degen- 
erate linked-diagram perturbation theory [ill, starting 
from the N'^LOW nucleon-nucleon potential [J]. The lat- 
ter is a low-momentum potential derived from chiral per- 
turbation theory at next-to-next-to-next-to-leading order 
with a sharp momentum cutoff at 2.1 fm~^. More ex- 
plicitly, we have derived H^s using the well-known Q- 



box plus folded-diagram method |ll| . where the Q-box 
is a collection of irreducible valence-linked Goldstone di- 
agrams, which we have calculated through third order in 
the N^LOW potential. 



TABLE I: Theoretical shell-model single-particle energies (in 
MeV) employed in present work (see text for details). 



nlj 



Single-particle energies 



Od| 
Od. 



■5/2 
■3/2 



lSl/2 



0.601 
5.121 
-0.793 



The hamiltonian Hcs contains one-body contributions, 
whose collection is the so-called S'-box [l2|. In realistic 
shell-model calculations it is customary to use a subtrac- 
tion procedure, so that only the two-body terms of HoS, 
which make up the effective interaction T^ff , are retained 
while the single-particle energies are taken from experi- 
ment. In this work, however, we have adopted a differ- 
ent approach employing single-particle energies obtained 



from the S'-box calculation. In this connection, it is worth 
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In Table 



noting that the observed lowest J" 
in ^^C are not pure single-particle states 
U our calculated single-particle energies are reported. It 
should be noted that the experimental one-neutron sep- 
aration energy of ^^C (S'„=1.22 MeV [Hi) is slightly un- 
derestimated. For the sake of completeness, in Table ?? 
we also report the two-body matrix elements of Vcs- 



TABLE II: Proton-proton, neutron-neutron, and 
proton-neutron matrix elements (in MeV). They 
are antisymmetrized, and normalized by a factor 
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We have performed calculations for carbon isotopes 
with A ranging from 16 up to 24, i.e. for systems with 
valence neutrons from A'vai=2 up to 10. In Fig. [l]the cal- 
culated ground-state energies of even- mass isotopes (con- 
tinuous black line) relative to ^^C are com par ed with the 
experimental ones (continuous red line) [l4| . The ex- 



perimental behavior is well reproduced, in particular our 
results confirm that ^^C is the last bound isotope; its 
calculated S2n is 601 keV to be compared with the eval- 
uation of 420 keV [IJ]. Moreover, our calculations predict 
that ^^C is unstable against one- neutron decay, the the- 
oretical Sn being -1.6 MeV. Therefore, our results fit the 
picture of ^^C as a Borromean nucleus. 
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FIG. 1: (Color online) Experimental [J, [iJI and calculated 
ground-state energies for carbon isotopes from A = 16 to 24. 
TV is the number of neutrons. See text for details. 



From the inspection of Fig. [U it can be seen that our 
calculations underestimate the experimental data. It is 
worth noting that this discrepancy may be "healed" by 
downshifting the single-particle spectrum so as to repro- 
duce the experimental ground-state energy of ^^C relative 
to ^^C. The results obtained with this downshift (—427 
keV) are reported in Fig. [l]by the black dashed line. 

In Fig. [5J we report the experimental and calculated 
low-energy levels of the odd-mass nuclei ^^C and ^^C. It 
should be pointed out that the experimental levels shown 
in Fig. [2]are the only observed bound states. We see that 
the theory provides a satisfactory description of these 

5 + 
2 



states, apart from the inversion between the J'" 

and ^ states in ^'^C. 

We have calculated the magnetic dipole moments of 
the J"" = (i+)i, (|+)i states in ^^c, and J"" = (|"^)i 
state in ^^C using an effective operator obtained at 
third order in perturbation theory, consistently with the 



derivation of Hos- We have obtained /z( 
— 1.803 n.m., and /i( 



n.m., ^(U 



21 



21 



-1.920 

-0.807 



n.m., to be compared with the experimental values 
±1.720 ± 0.009 [li n.m., ±1.758 ± 0.030 n.m. O, and 
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FIG. 2: Experimental [lg| and calculated low-energy spectra 
for ^''C and ^^C. 



Two magnetic dipole transition rates have been mea- 
sured in "C: the B{M1; \^ -^ % 



and B(Af 1; |+ 



I ). Our calculated values are respectively 0.0941 ^^ 
and 0.095 //^, to be compared with the experimental val- 
ues of (0.0010 ±0.0001) ix\ and (0.082 t^o^s) M?,r [Hi- 
lt is worth mentioning that the anomalous quenching of 
the strength observed in the \ ~^ | transition has 
been recently reproduced in Ref. [7|] in the framework 
of a shell-model calculation modifying the tensor com- 
ponents and the T = 1 monopole terms of the empirical 
effective Hamiltonian SFO j20|| defined in the the full -psd 
model space. 



TABLE III: Effective reduced single-neutron matrix elements 
of the magnetic dipole operator Ml (in n.m.). 



nJaja nbhjb 



(a||Afl||b) 



Odg/a Od5/2 

04/2 0^3/2 
Odi/2 0^5/2 
0^3/2 Od3/2 
0^3/2 lSl/2 
lSi/2 0^3/2 
lSl/2 lSi/2 



-2.553 
2.743 
-2.743 
1.414 
0.179 
-0.179 
-2.298 



In Tables IIIII we report the effective reduced single- 
neutron matrix elements of the Ml operator. 

In Fig. [3] we report the experimental excitation ener- 
gies of the yrast 2+ states as a function of A and compare 
them with our calculated values. It can be seen that the 
observed energies are reproduced nicely. We also report 
our predicted excitation energy, 4.661 MeV, for the un- 
bound ,r = 2^ state in ^^C. 



TABLE IV: Effective reduced single-neutron matrix elements 
of the electric quadrupole operator E2 (in c fm^). 
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lSl/2 
04/2 
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-3.357 
-1.866 
-1.826 
1.866 
-2.803 
1.434 
-1.826 
-1.434 



5.6 








4.8 


• Expt. 

• Calc. 




. 


4.0 






/ 


3.2 






/ 


2.4 






^ 


1.6 


•— — __ 


* 









10 



12 



14 



16 



FIG. 3: (Color online) Experimental [2l|, [23] and calculated 
excitation energies of the yrast J^ = 2^ states for carbon 
isotopes from A = 16 to 22. A'' is the number of neutrons. 



The behavior of the theoretical 2+ excitation energies 
confirms that there is no A^ = 14 subshell closure. We 
obtain a ground-state wave function in ^"^C with only 
14% of the (1^^5/2)^ configuration (calculated using the 
OXBASH shell- model code |2j|). This is a direct conse- 
quence of the fact that the effective single-particle energy 
|24| of the lsi/2 state is the lowest one all along the car- 
bon isotopic chain, as can be seen in Fig. |31 One can 
infer that the large energy gap between the 0^3/2 level 
and the 0^5/2, lsi/2 levels is responsible for the iV = 16 
subshell closure. 

It is worth mentioning that the behavior of the yrast 
2^ state in C isotopes is quite different from that shown 
by O isotopes, where a TV = 14 subshell closure has been 
evidenced. This is clearly related to the removal of the 
two protons from the Opi /2 level, which, in the case of 
oxygen isotopes, interact with the sd neutrons giving rise 
to a downshift of the 0^5/2 neutron level relative to the 
lsi/2 one |8!|. In this connection, we have verified that, 
using our effective hamiltonian, the effective single par- 
ticle energy of the 0^5/2 state in ^^O is lowered by 2.1 
MeV relative to that of the lsi/2 one. 




FIG. 4: (Color online) EfTective single-particle energies of the 
neutron lsi/2, Od^/2, and 0^3/2 orbits from ^ = 16 to 22. A'' 
is the number of neutrons. 



We have calculated the B{E2;2^ — > Of) transition 
rates up to ^°C employing the effective operator of Ta- 
ble IfVI Our results are compared with the experimental 
data in Table |Vl We see that the agreement is quite 
good, providing evidence for the reliability of our calcu- 
lated effective operator which takes into account micro- 
scopically core-polarization effects up to third order in 
the NN potential. It is well known [lOJ, l2l|, 12^ |25| that 
these observed transition strengths are strongly hindered 
when compared with the values obtained by an empiri- 
cal formula based on the liquid-drop model [26|. In this 
connection, we would like to point out that our effective 
reduced single-neutron matrix elements (see Table lIVp 



correspond to a neutron effective charge of about 0.4e 
with an harmonic-oscillator parameter b — 1.72 fm. 



TABLE V: Experimental and calculated reduced transition 
probabilities B{E2\ 2t — ^ 0^^) (ine^fm*). Superscripts a, h, c 
refer to [ly| , [23| , and [22| , respectively. The experimental er- 
rors reported are the statistical and systematic ones, respec- 
tively. 



Nucleus 



Calc. 



Expt. 



^C 



^C 



1.8 


2.6 ±0.2 ±0.7 




4.15 ±0.73 ^ 


3.0 


4.3 ±0.2 ±1.0 


3.7 


<3.7" 



In summary, we have given here a shell-model descrip- 
tion of heavy carbon isotopes, using a fully microscopic 
approach. This has been done by deriving from the realis- 
tic chiral NN potential N'^LOW both the single-particle 
energies and residual two-body interaction of the effective 
shell-model hamiltonian. Our calculations have led to a 
sound description of these isotopes when approaching the 
neutron dripline. As a matter of fact, they confirm the 
disappearance of the TV = 14 subshell closure that is 
present in the oxygen chain, and predict the iV = 16 one. 
In particular, we reproduce successfully the fact that ^^C 
is the last bound isotope. It is worth mentioning that 
this nucleus, which has only recently been observed [ij, 
is one of the most exotic of the 3000 known isotopes, its 
N/Z ratio being 2.67. The quality of the results obtained 
makes us confident that our approach may be a valuable 
tool to study nuclei far from the valley of stability. 
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